Measurement of the Xci Polarization in il){2S) — > 7Xc2 
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The polarization of the Xc2 produced in 4>{2S) decays into 7Xc2 is measured using a sample of 
14 X lO'^ ^{23) events collected by BESII at the BEPC. A fit to the Xc2 production and decay angular 
distributions in 'tp{2S) 7Xc2, Xc2 tt+tt" and K^K~ yields values x = A\/Ao — 2.08±0.44 and 
y — A2/A0 — 3.03±0.66, with a correlation p = 0.92 between them, where ylo,i,2 are the Xc2 helicity 
amplitudes. The measurement agrees with a pure El transition, and M2 and E3 contributions do 
not differ significantly from zero. 

PACS numbers: 13.20.Gd, 13.25.Gv, 13.40.Hq, 14.40.Gx 

I. INTRODUCTION II II H- In general, it is believed that V (25") -> 7Xcj 



The radiative transition between charmonium states 
has been studied extensively by many authors Q, |^ |^ 
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is dominated by the El transition, but with some M2 
(for Xci and Xc2) and E3 (for Xc2) contributions due to 
the relativistic correction. These contributions have been 
used to explain the big differences between the calculated 
pure El transition rates and the experimental results 
They will also affect the angular distribution of the ra- 
diative photon. Thus the measurement of the angular 
distribution may be used to determine the contributions 
of the higher multipoles in the transition. 

Furthermore, for ip{2S) — > 7Xc2, the E3 amplitude is 
directly connected with Z?-state mixing in ipi'^S) which 
has been regarded as a possible explanation of the large 
leptonic annihilation rate of -0(3770) 0. Since recent 
studies 0, 13 also suggest the S- and D-wave mixing of 
'0(25') and ■0(3770) may be the key to solve the long- 
standing "pTT puzzle" and to explain 0j(377O) non-DD 
decays, the experimental information on multipole am- 
plitudes gains renewed interest. 

Decay angular distributions were studied in 0(25) — > 
7Xc2 by^ the Crystal Ball experiment using 0(25) 
j'^J/ip [of; the contribution of the higher multipoles were 
not found to be significant but the errors were large 
due to the limited statistics. In the present analysis, 
0(25) 7Xc2 — ^ 77r+7r~ or jK~^K~ decays will be used 
for a similar study. The analysis on these channels has 
the advantage that there is no background from Xci since 
the Xci ~^ TT^TT" and K~ processes are forbidden by 
parity conservation. 



II. THE BES EXPERIMENT 

The data used for this analysis are taken with the 
BESII detector at the BEPC storage ring operating at 
the 0(25;V The number of 0(25) events is 14.0 ± 0.6 
million [lOj, determined from the number of inclusive 
hadrons. 

The Beijing Spectrometer (BES) detector is a con- 
ventional solenoidal magnet detector that is described 
in detail in Ref. [Till- BESII is the upgraded version of 
the BES detector |ia|. A 12-layer vertex chamber (VC) 
surrounding the beam pipe provides trigger information. 
A forty-layer main drift chamber (MDC), located ra- 
dially outside the VC, provides trajectory and energy 
loss (dE/dx) information for charged tracks over 85% 
of the total solid angle. The momentum resolution is 
ap/p = 0.017-\/l +p'^ {p in GeV/c), and the dE/dx reso- 
lution for hadron tracks is ^ 8%. An array of 48 scintilla- 
tion comitcrs surrounding the MDC measures the time- 
of-flight (TOP) of charged tracks with a resolution of 
^ 200 ps for hadrons. Radially outside the TOP system 
is a 12 radiation length, lead-gas barrel shower counter 
(BSC). This measures the energies of electrons and pho- 
tons over ~ 80% of the total solid angle with an energy 
resolution of as/E = 22%/VE {E in GeV). Outside of 
the solenoidal coil, which provides a 0.4 Tesla magnetic 
field over the tracking volume, is an iron flux return that 
is instrumented with three double layers of counters that 



identify muons of momentum greater than 0.5 GeV/c. 

A GEANT3 based Monte Carlo (MC) program with 
detailed consideration of detector performance (such as 
dead electronic channels) is used to simulate the BESII 
detector. The consistency between data and Monte Carlo 
has been carefully checked in many high purity physics 
channels, and the agreement is quite reasonable. 

MC samples of 0(25) 7XcO,2 77r+7r^ and 
0(25) 7XcO,2 jK'^K~ are generated according to 
phase space to determine normalization factors in the 
partial wave analysis. MC samples of e+e^ — s- (7)e^e~, 
^^(25) ^ (7)e+e-, e+e" ^ (7)/^+^", 0(2^) -> 
(7)/i+M-, and_0(25) ^ XJ/^, J/0 ^ (7)/^+^" {X ^ 
77, TT^TT^, 7r+7r~, and ry) are used for background estima- 
tion. 



III. EVENT SELECTION 

For the decay channels of interest, there are two high 
momentum charged tracks and one low energy photon. 
The candidate events are required to satisfy the following 
selection criteria: 

1. At least one photon candidate is required. A neu- 
tral cluster is considered to be a photon candidate 
when the angle between the nearest charged track 
and the cluster in the xy plane is greater than 15°, 
the first hit is in the beginning six radiation lengths 
of the BSC, and the angle between the cluster de- 
velopment direction in the BSC and the photon 
emission direction in the xy plane is less than 37°. 
There is no restriction on the number of extra pho- 
tons. 

2. Two good charged tracks with net charge zero are 
required. Both tracks must satisfy |cos0| < 0.65 , 
where 9 is the polar angle of the track in the lab- 
oratory system. This angular region allows use of 
the fi counter information to eliminate p,~^ n~ back- 
ground. 

3. To remove Bhabha events, the total energy de- 
posited in the BSC energy by the two charged 
tracks is required to be less than 1 GeV, or 

dE/dx (dE/dx),„aas~{dE/dx)^^„ r i, x i • 

Xe — - — - — ^ — - — for each track is 
required to be less than -3. Here {dE/dx)nicas and 
{dE / dx)c^p are the measured and expected dE/dx 
energy losses for electrons, respectively, and a is 
the experimental dE/dx resolution. This removes 
almost all events with two electron tracks but keeps 
the efficiency high for the signal channels. 

4. To remove /i^/i^ backgrounds, MUID^ + 
MUID^ < 3 in the 77r+7r- channel and MUID+ + 
MUID^ < 5 in the -fK+K~ channel are re- 
quired. Here MUID is the number of ji counter 
hits matched with the MDC track and ranges from 
to 3. "0" means not a /i track, while "1", "2" 
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and "3" means a loose, medium, or strong /i can- 
didate [3. 

5. To remove cosmic rays, \t^oF ~ ^tof\ < 4 ns is 
required, where Itof is the time recorded by the 
TOF. This removes all the cosmic ray events with 
almost 100% efficiency for the channels of interest. 

6. Four-constraint kinematic fits are performed with 
the two charged tracks and the photon candidate 
with the largest BSC energy under the hypothe- 
ses that the two charged tracks are either tt+tt" or 
K'^K^ , and the kinematic chisquares, xi E^^d Xk, 
are determined. If < Xk ^^'^ confidence 
level of the fit to •4'{2S) 77r"'"7r~ is greater than 
1%, the event is categorized as 77r+7r^; otherwise, 
if Xk < xi Etnd the confidence level of the fit to 
tp{2S) ^K~^K~ is greater than 1%, the event is 
categorized as jK^K~. 

After imposing the above requirements, the invariant 
mass distributions for the selected 77r+7r~ and ^K^K^ 
candidates are shown in Fig.^ Clear Xco and Xc2 signals 
can be seen while the background level is low. 

Simulated background events passing the selection cri- 
teria for the 77r+7r" and ^K'^K" channels are also plot- 
ted in Fig.Q] The excess background in the 77r+7r^ mode 
near 3.7 GeV/c^ is due to the large ip{2S) tt+tt" 
branching ratio from the PDG The backgrounds 

under the signal regions are (7) and e+e~ (7) events 

either from QED processes or from ^{2S) decays. 

Requiring the invariant mass of the two charged tracks 
be between 3.54 and 3.57 GeV/c^ to select Xc2, 418 
77r"'"7r~ events and 303 jK~^K~ events are selected. The 
fractions of background are (1.6 ± 0.5)% for 77r+7r^ and 
(2.8±0.6)% for jK+K-, as estimated from Monte Carlo 
simulation, in agreement with the expectation from the 
measured misidentification efficiencies in data. 

Monte Carlo simulation also determines that the 
K~ contamination in the tt+tt^ sample is about 9%, 
and the tt^tt^ contamination in the K^K~ sample is 
about 34%. The effect of the cross contamination on the 
fit of the helicity amplitudes will be discussed later. 

IV. THE FIT OF THE HELICITY AMPLITUDES 

The ip(2S) — > 7Xc2 helicity amplitudes are determined 
by a maximum likelihood fit to the decay angular distri- 
bution [Hill 

W2{e-,,eM, (t>M) = ix^ sin^ 9^ sin^ 26*4^ + (1) 
(1 + cos^ e^)[{?, cos^ Om - if + sin'' Om] + 
Vix sin 20^ sin 26* j\/ [3 cos^ Om — 1 — 
-%/6?/ sin^ 9m] cos 0a/ 
+V&y SIT? 6^ sir? 6m (3 cos^ 6m — 1) cos 20^/ , 



where x = Ai/Aq, y = A2/AQ, ^0,1,2 are the Xc2 helicity 
amplitudes, 6^ is the polar angle of the photon in the 
laboratory system, and 6m and (pM are the polar and 
azimuthal angles of one of the mesons in the Xc2 rest 
frame with respect to the 7 direction. (pM = is defined 
by the electron beam direction. 

Fitting the 77r+7r^ and ^K^K^ data, we obtain 

x^ = 1.97 ±0.64, = 3.03 ± 1.07, = 0.96, 
XK = 1.77 ±0.54, = 2.36 ±0.82, p/f = 0.94, 

where the errors are statistical and p^rj Pk are the corre- 
lation factors between x and y for 77r+7r^ and jK^K^, 
respectively. The comparison between data and the fit 
is shown in Fig. |21 Good agreement is observed in all 
angular distributions for both the "/tt^tt^ and "fK^ 
channels. 

Since the value of the likelihood function does not pro- 
vide a measurement of the goodness of fit, Pearson's x^ 
test is used. The data are divided into 3 x 3 x 4 = 36 bins 
in cosOj, cos 6 m and </>m- The x^ is calculated using 

(„DT „MC\2 

where nf^ is the observed number of events in the ith bin 
and nf^ is the corresponding number of events predicted 
by Monte Carlo using x and y fixed to the values deter- 
mined in this analysis. We obtain x^/ndf — 30.19/35 = 
0.86 and xji/ndf = 43.57/35 = 1.24 for the 77r+7r- and 
^K~^K~ channels, respectively, where ndf is the number 
of the degree of freedom. These results show that the fits 
are good. 

V. ERROR ANALYSIS 
A. Input output checking 

The fitting procedure is tested using Monte Carlo simu- 
lated samples. With input parameters Xm — VS ~ 1.732 
and yin = VG ~ 2.449, fitting a Monte Carlo sample of 
50,000 selected events gives the results Xout — 1.74±0.04, 
yout = 2.45 ±0.07, and p = 0.94, which are in good agree- 
ment with the input values, indicating the validity of the 
fitting procedure. 

Dividing the 50,000 events into 100 subsets of 500 
events each (about the same size as the real data sample) 
gives the distribution of fitting results shown in Fig. 13 x 
and y are positively correlated, and the fitting results are 
distributed in a relatively broad area due to the limited 
statistics of the subsets. 



B. Systematic errors 

Systematic errors from background, from the 77r+7r^ 
and ^K^K~ cross contamination, from the Monte Carlo 
simulation of the detector response, etc. are considered. 
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FIG. 1: Invariant mass distributions of the two charged tracks in (a) jtt^-k and (b) 'yK'^K . Dots with error bars are data, 
and the shaded histograms are the MC simulated backgrounds. 




FIG. 2: Comparison between data and the final fit for 
77r"'"7r~(left) and yK^K~ (right), where dots with error bars 
are data and the histograms are the fit. 



1. Background contamination 

Backgrounds remaining after event selection are 
11^ jjT {'y) and e^e~(7) events, and the fractions of back- 
grounds in 77r+7r~ and 'yK~^K~ channels are estimated 
by Monte Carlo simulation and checked with data. In 
the fit, background is not considered, but the effect on 
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x=1 .74±0.04 
y=2.45±0.07 



FIG. 3: Distribution of fitting results for x and y for Monte 
Carlo simulated samples. The black dot with error bar is the 
result for all 50,000 events. The circles are the fitting results 
for the subsets after dividing the sample into 100 subsets. 



the helicity amplitudes is estimated using Monte Carlo 
simulation. By adding the amount of MC background 
mentioned in Sec. IIIII into the pure MC sample, the fit 
yields shifts of the fit results. These shifts are taken as 
corrections to the results obtained from data. By varying 
the background fraction in the fit, the uncertainty due to 
the background contamination can also be determined. 
It is found that the corrections to the 77r+7r~ results are 



0.19 ±0.04, A 



Try 



Ajf, = 0.25 ± 0.04, AKy 



0.34 ± 0.07; and for jK+R- 
= 0.47 ±0.08. 
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2. "^■k'^-k and jK'^K cross contamination 

In order to study the error from 77r+7r^ and ^K^K^ 
cross contamination, Monte Carlo samples of jtt^tt^ and 
jK~^K~ with X = -\/3, y = \/6 are generated and mixed 
according to the amount of cross contamination deter- 
mined by Monte Carlo simulation. It is found that the re- 
sults from this mixed sample are mostly unchanged from 
those of the pure Monte Carlo sample, even when the con- 
tamination is doubled. This is understandable since the 
angular distributions of jTr^n^ and jK^ K~ are identi- 
cal. From the comparisons of many Monte Carlo samples 
with different fractions of cross contamination, the errors 
on X and y are determined to be 0.01 and 0.06 for Xj^ and 
j/tt, and 0.10 and 0.14 for xk and yx, respectively. 

3. MC simulation of the detector response 

The consistency between data and the Monte Carlo 
simulation of the detector response for Xc2 events can be 
determined using Xco events, although the absolute angu- 
lar distributions are different. The angular distribution 
of XcO decays is unambiguous, i.e. Wo = 1-l-cos^ 9^. Note 
that Eq. 1 with the (3cos^6'm — 1) term replaced by 1 
is equal to Wq when both x and y are zero. Therefore, 
if we fit the angular distribution of XcO the same as Xc2 
using the modified Eq. ^ x and y should be 0. The dif- 
ference from zero gives the systematic error due to the 
MC simulation of the detector response. For the 77r+7r~ 
channel, 0.18, 0.05, and 0.24 are obtained for x^^, y^ and 
Ptt, respectively, and for the jK~^K~ channel, 0.13, 0.08, 
and -0.24 are obtained for xk-, yK and pK- The results 
are dominated by the statistical errors of the fit due to 
the limited XcO samples, although they are already much 
larger than the corresponding Xc2 samples. Comparisons 
between data and Monte Carlo simulation for XcO events 
are shown in Fig. ^ good agreement is observed. 

4. Other sources 

Other sources of error are from systematic errors asso- 
ciated with the simulation of the mass resolution of the 
Xc2, the photon detection efficiency, the MDC tracking ef- 
ficiency, the kinematic fit, the total number of the '0(25') 
events, the trigger efficiency, etc. These systematic er- 
rors will affect a branching ratio measurement, but will 
not affect the measurement of the angular distribution. 
Their effects on the helicity amplitude measurements are 
neglected. 

5. Total systematic error 

The systematic errors and the correlation factors from 
all the above sources are listed in Table Here the cor- 
relation factors (p^ and pk) from background contami- 




FIG. 4: Comparison of angular distributions between data 
(dots with error bars) and Monte Carlo simulation (his- 
tograms) for XcO tt^tt" (left) and XcO — > K'^ K~ (right). 
Fitting the XcO angular distributions provides a way to esti- 
mate the systematic error due to the Monte Carlo simulation 
of the detector response. 



nation and 77r+7r^ and jK^ cross contamination are 
set to 1, and the total correlation factor p is calculated 
with p = 5]^£i£kL2jii^ where i runs over all the systematic 

errors. The total systematic errors are 0.19 and 0.11 for 
X and y in 77r+7r^ and 0.17 and 0.18 in 'yK^K^ . 



TABLE I: Summary of the systematic errors and correlations. 



Source 


^TT p7V 


Xk Vk Pk 


Background contamination 
■K / K cross contamination 
MC simulation 


0.04 0.07 1 
0.01 0.06 1 
0.18 0.05 0.24 


0.04 0.08 1 
0.10 0.14 1 
0.13 0.08 -0.24 


Total 


0.19 0.11 0.29 


0.17 0.18 0.49 



VI. RESULTS AND DISCUSSION 

After applying the corrections due to the background 
contamination, we obtain 

x^ = 2.16 ±0.64 ±0.19, 
= 3.37±1.07±0.11, 
°* = 0.96, p'y' = 0.29 
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from 77r"'"7r and 

XK = 2.02 ±0.54 ±0.17, 
Vk = 2.83 ±0.82 ±0.18, 

from ^K'^K~ , where the first errors are statistical and 
the second are systematic, and and p^^^^ are the 

correlation factors between x and y of the statistical and 
systematic errors. 

Combining the statistical and systematic errors yields: 

x^ ^ 2.16 ±0.67, = 3.37 ± 1.08, = 0.93, 
XK = 2.02 ±0.57, = 2.83 ±0.84, = 0.91. 

The results from 77r"'"7r~ and jK^K^ are in good agree- 
ment. Combining them, we obtain 

X = 2.08 ±0.44, 
y = 3.03 ±0.66, 
p ^ 0.92. 

The combination assumes no correlation between 77r+7r^ 
and jK~^K~ for both statistical and systematic errors. 

Comparing with the measurement obtained by the 
Crystal Ball 9], 4(V/ ^ 7X) = 0.132t°:°?«, 
this measurement gives the quadrapole amplitude 
a2 = — 0.051!too36 and the octupole amplitude = 
-0.027lJJ;!^29 El- Neither result significantly differs from 
zero. The results are in good agreement with what is 
expected for a pure El transition. As for the D-state 



mixing of ■0(25'), our results do not contradict the pre- 
vious theoretical calculation within one standard devia- 
tion [13. 



VII. SUMMARY 

The helicity amplitudes of 'ip{2S) — > 7Xc2 are measured 
for Xc2 TT+TT" and K+R-, and x = 2.08 ± 0.44, y = 
3.03 ± 0.66 with correlation p = 0.92 are obtained. The 
results are in good agreement with a pure El transition, 
but still do not have the precision to strongly limit the 
higher multipoles. 
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